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Abstract 

The A3(L) statistic of Random Matrix Theory is defined as the average of a set of random num- 
bers {5}, derived from a spectrum. The distribution p{6) of these random numbers is used as the 
basis of a maximum hkehhood method to gauge the fraction x of levels missed in an experimental 
spectrum. The method is tested on an ensemble of depleted spectra from the gaussian orthogonal 
ensemble (GOE) , and accurately returned the correct fraction of missed levels. Neutron resonance 
data and acoustic spectra of an aluminum block were analyzed. All results were compared with an 
analysis based on an established expression for A-^{L) for a depleted GOE spectrum. The effects 
of intruder levels is examined, and seen to be very similar to that of missed levels. Shell model 
spectra were seen to give the same p{S) as the GOE. 
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I. INTRODUCTION 



Neutron resonance data provide us with a high resolution picture of the eigenvalues of the 
nuclear hamiltonian at high excitation energies. This was the birthplace and testing ground 
for Random Matrix Theory (RMT) as a model for quantum chaos. For a brief history of 
RMT see [l], and for a review of RMT and nuclear structure, see The correspondence 



between the fluctuation properties of nuclear spectra and those o 

n 

nsemble has been verified many times in neutron resonances p 



the Gaussian Orthogonal 
and proton resonances 



3]. Furthermore, shell model calculations exhibit many of the fluctuation properties of the 
GOE js, 9|. For an account of tests of RMT in nuclear physics see jlO |. 

The question of the completeness of an experimental spectrum is important. One needs 
a gauge of the fraction, x, of the levels missed in a given experimental spectrum. RMT 
has already been used to this end. The fraction of levels not observed due to the finite 
resolution and sensitivity of the detectors will chan ge t he distribution of widths from the 
Porter Thomas distribution which follows from RMT ll| . The nearest neighbor distribution 



(nnd) is another commonly used statistic. The nnd for a pure spectrum follows the Wigner 
distribution, 



Pis) 



-se 



-TTS^/4: 



(1) 



where s = S/D, S being the spacing between adjacent levels, and D is the average spacing. 
The nnd of a spectrum incomplete by a fraction x is given by 



P(s) = Y.{l~x)x^P{k;s) 

k=0 



(2) 



where P{k; s) is the k nearest neighbor spacing, -Efc+j — Ei. This was first introduced as 
an ansatz in j?], and rederived in |12|] and 13 1. Eq. |2]was used by Agvaanluvsan et al as 
ihe basis for a maximum likelihood method (MLM) to determine x for incomplete spectra 
12|. The A3(L) statistic (also called the spectral rigidity) introduced by Dyson [ij] is a 
commonly used statistic. It is defined as a spectral average: 



A3(L) = /mmA,B j J^^^"^ dE' [ M{E') - AE' - Bf 



mL)) 



(3) 



where J\f{E) is the cumulative level number, the number of levels with energy < E' ( its slope 
is the level density p(-E) ). A and B are chosen to minimize Sl{L). They are recalculated for 



each i. A series of evenly spaced levels would make M{E) a regular staircase, and ^^{L) = 
Y^. At the other extreme, a classically regular system will lead to a quantum mechanical 
spectrum with no level repulsion, the fluctuations will be far greater, and A^^L) = ^. The 
angle brackets mean the average is to be taken over all positions Ei of the window of length 
L. 

An analysis of neutron resonance data using the A3(L) statistic and the MLM of Ag- 



vaanluvsan et al was performed in [15| and 



16| with consistent results. When both methods 



were tested on ensembles of depleted GOE spectra, the mean values of x were correct, but 
the uncertainties were large, for realistic spectrum sizes. 

In this paper we present a new method to test experimental data for missed levels. It is a 
MLM based on the definition of the A3(L) statistic. Instead of concentrating on the spectral 
average of the random numbers in Eq. [3l we consider instead their distribution. If 

there are D levels in the spectrum, then, allowing for setting the zero of the energy scale at 
the lowest level, there are D — L — 1 values of Sl{L) for each value of L. This amounts to a 
large sample size of random numbers. In this paper we use numerical simulation to get an 
expression for the distribution of Sl{L) for depleted spectra with the fraction missed, x, as 
a parameter, and base a MLM on this distribution. The method will return a most likely 
value of X for each L. 

In 13| Bohigas and Pato gave an expression is given for A3(L) statistic for incomplete 



spectra. The fraction of missed levels x is both a scaling factor and a weighting factor and 
A2{L,x) is the sum of the GOE and poissonian result: 



, ,^ x2aG0E, 



A3(L, x)=x'^ + {l- xYAY^'^iL/il - x)). (4) 

The A3(L) statistic of an experimental spectrum can be compared with this expression and 
the best x found. We will see however, that the uncertainty in x is large for this method. 

In the next section we describe the process of making, unfolding and depleting GOE 
spectra, and the calculation of Sl{L). In Sec. [1111 we discuss the cumulative distribution 
function, Af{6), of Sl{L) for spectra with depletion x. A three parameter fit was sufficient 
for each x. The parameters were found for 0.00 < x < 0.30 in steps of 0.01. We then fit the 
parameters as functions of x. Now we have A/'(5) with x as a continuous parameter. This is 
the basis for our MLM. In Sec. [IV]the MLM is developed and tested on ensembles of depleted 
GOE spectra. In Sec. |V]we tested Eq. |4]and used it to return x for depleted GOE spectra, 
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and compared these results with the MLM. In Sec. IVII the method is apphed to neutron 
resonance data and ultrasonic spectra measured from an aluminum block. The results are 
compared with previous investigations. The effect of intruder levels on A3(L) and p(S) is 



very similar to missed levels. This issue is addressed in Sec. IVIII Recent developments 17 1 
questioned the validity of using RMT to model fluctuations of complex spectra. To see if 
N'iS) could discriminate between the GOE and more physical models we looked, in Sec JVIIIi 
at shell model spectra. In Sec. IIXI we make some concluding remarks. 

II. RMT CALCULATIONS 

To do our RMT analysis we need to generate an ensemble of random matrices, diagonalize 
them, unfold them, deplete them, and calculate Sl{L) for each of them. We will give a brief 



description of this process here, it is described in morbid detail in |15|. 

The appropriate ensemble for this analysis is the Gaussian Orthogonal Ensemble (GOE) 
as it describes real, time-reversal-invariant systems. This ensemble is the set of random 
matrices H whose elements are normally distributed matrix elements, Hij, having 

for the off-diagonal and diagonal elements respectively. The width of the distribution is 
arbitrary, and we choose a = 1. Each of these matrices is diagonalized. We made an 
ensemble of 3000 matrices, each of dimension D = 3000. Each matrix has an approximately 
semicircular level density, with p{E) = ^JAN — E^, for l-E] < 2\/N and p{E) = otherwise. 
Interesting as this may be, it is not germane to our analysis. The currency of RMT is 
fluctuations, and in order to compare GOE results with experimental data we must remove 
this long range (secular) structure from all spectra, including the experimental data, with a 
process called unfolding. The basic idea is to rescale the energy axis to give a uniform level 
density of one level per unit energy, on average. To work with smaller spectra we just take 
a section from the large {D = 3000) spectra of whatever size we want. 

In a spectrum with picket fence of levels, spaced 1 unit apart, like the harmonic oscillator 
spectrum, Af{E) is a staircase with steps 1 unit high and 1 unit long, and the spectrum is 
said to be rigid. An arbitrary spectrum has N'{E) = i, Ei < E < -Ej+i. The A3(L) statistic 
is a measure of fluctuations of Af{E) from a regular staircase, and its definition is the square 
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of the difference between this stairs and a straight hne. In the harmonic oscillator case Eq. 
[3] can be integrated directly to get A3(L) = 1/12. The situation is messier for an arbitrary 
spectrum. Using Af{E) = i, Ei < E < -Ej+i, in Eq. (E]), and performing the integral between 
two adjacent levels, we come to 

I i+L-l i-Ej+l 



^W = - Y: I ^ dE'U-AE'-Bf 
L 7~ JEi 



= -x {C + VA^ + WA + XAB + YB + ZB"^), 
L 

where C = E'jtt' f{E,^i ' E,), V = \{El^ - Ef), W = ^fi-' -j{E]_,, - Ej), X = 
{Ef^L - Ef), Y = -2j{Ej+i - Ej), Z = {Ei+L - Ei). A and B need to minimize 

5\{L), and this leads to the constraints d{/\\)ldA = and d{Al)/dB = 0. These equations 
readily give A and B that minimize A3(L) as follows: 

XF - 2WZ _ WX - 2VY 
4VZ-X^ ' 4VZ-X^ ' 

To generate spectra with specific values of D and x, we take a section of levels 
from the middle of an unfolded GOE spectra, randomly drop a fraction x of them, and 
then contract this spectrum by a factor of 1 — x to restore a level density of 1 level per 
unit energy. So this gives D levels "detected" from the original spectrum, and a fraction 
X "missed". Note that x is not a continuous parameter. Consider an experimental run of 
124 levels. The true spectrum could have 125, 126 . . .or 130 levels, there being 1, 2, . . .or 
6 levels missed, in which case x would have values of 0.8% , 1.59% , 2.36% , 3.13% , 3.88%, 
or 4.62%. This should be kept in mind when testing various schemes for estimating x. 

Now we can make ensembles with a range of x and D. We chose to restrict our ensembles 
to 1500 elements each, and chose D = 1000 and 0.00 < x < 0.30 in steps of 0.01. 



III. THE DISTRIBUTION OF 6i{L) 

Given an unfolded GOE spectrum of size D, where a fraction x has been missed, A3(L) 
will be the (spectral) average of the set of D — L — 1 numbers Sl{L). We would like to get 
the probability density of these numbers as a function of s. In what follows we will write 6 
for Sl{L), dropping all subscripts assuming as fixed value for L. So the probability density 
of Sl{L) will be written p{S), and the cumulative distribution will be written Af{6). 
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FIG. 1: (color online) The cumulative distribution function, J\f{5). Each set of 3 lines corresponds 
to, from left to right, x = 0.00, 0.04, and 0.08. The leftmost set of 3 lines of lines are for for L = 10, 
and the rightmost set have L = 40. The 3-parameter fit, EqlS]is plotted in black for each case, and 
is indistinguishable from the ensemble data (blue). 

We should explicitly state that this analysis assumes that p{6) is ergodic, in that the 
distribution of 6 is the same for one huge matrix as it is from the superposition of many 
small ones. It is true that the level spacing distribution p{s) is ergodic. A histogram of 
the 3000 spacings from a single D > 3001 spectrum is the same as a histogram of the 
s = EiQ — spacing of 3000 spectra with D > 16, for example. But A3(L) is defined as a 
spectral average, so while it seems obvious that it is ergodic, it seems prudent to state the 
assumption explicitly, given that that A3(L) can be very different from one spectrum to the 



next see Ref. 



15|. 



We proceeded by guessing the functional form of p{S), and were surprised to see that 
p{6) was a simple function of log 6 for a wide range of L and x, and this motivated us to 
parameterize the cumulative distribution function. We used following parameterization: 

Af{5) = ^(l-Erf[a + Mog5 + c(log5)2]). (5) 

This yields, on differentiation 

p{6) = ^exp[-(a + Mog5 + clog52) ]( + ^). 6) 



In Fig. [T] we see the ensemble average of J\f{6) for L = 10 and 40, with x = 0.00, 0.04 



and 0.08. The L = 40 case is more sensitive to depletion than the L = 10 case, one can 
clearly see that the spread in the x = 0.00, 0.04 and 0.08 lines is greater. This is reasonable 
because the bigger the window size, L, more likely it is to fall across the site of a missed 
level. Consider a spectrum with D=1000, and x = 0.02, there are 20 sites for missed levels. 
With L = 10, there are 989 positions for the window, and there are at most, 10 positions, i, 
where Sl{L) will be different from the x = 0.00 case, compared to 40 for the L = 40 case, so 
low L values give a less sensitive distribution. We note that there are just 3 parameters in 
the fits, shown in blue in the figure, even so their curves lie on the ensemble average values. 
The spread of the averaged values for each point in the J^{S) graphs was of order 10~^. 

The method used to numerically make the Af{S) is best illustrated by the following 
example for L = 20, and x = 0.03. Take a GOE spectra, with D = 3000. Take the middle 
100/0.97=1031 levels to avoid end effects. Unfold it. Randomly drop 31 levels. Contract 
the spectrum by a factor of (1 — x), now the level density is 1. Calculate the set of 1000-20- 
1=979 numbers {6}. Sort them. Do this 1200 times, and get {6} the average of the sorted 
sets. The standard deviation of the sets was ~ 10^'^. Now pair 5i, the i*'' element of {5}, 
with i/979, to get the set {6i,i/979}. A plot of {5i,i/979} is a graph of 7V((5). 

The values of a, b, and c were obtained by fitting Eq|5]to {6i, i}. This was done for L = 5 
to 90 in steps of 5, and for x =0.00 to 0.30 in steps of 0.01. The parameters were smooth 
values of x for all L. See FigjJl For each value of L the parameters were fit to smooth 
functions of x, aL{x) = ao + aiy/x + aix + a2x'^ , with similar expressions for biix) and cl{x). 
The values of oq, ai, oi, 02, ... , ci, C2 were calculated for all values of L from 5 to 100. We 
now have a probability density for S^^{L) with continuous parameter: 

p{6,x) = — ^exp[-(aL(x) + 6i(x) log 5 + cl{x) log 6^) ]{-^+ , )• (7) 

y/n 

In Figj3]we have some examples of p{5, x). When the fitted values, a{x), b{x), and c{x) were 
used in p{6, x) the results were indistinguishable from when those values of a, b, and c that 
were got from the fitting procedure were used. Again, we anticipate from the graph that 
higher values of L will be more useful in gauging x. In the next section we see how these 
parameters are used to get x for an unfolded spectrum. 

As a test of our machinery, we checked that p{6) was independent of x for completely 
uncorrelated (poissonian) spectra, and it was. 
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FIG. 2: The parameters o, b, and c of EqjSj On the left panels, the parameters are plotted vs x 
for fixed values of L. On the right panels we have plots of the parameters vs L for fixed x. From 
data like this we extracted aL{x), bL{x) and cl{x) 

IV. THE MAXIMUM LIKELIHOOD METHOD. 

Now that we have the probability distribution of 5\{L) parameterized, Eq. [TJ we can use 
it to find the most Hkely value of x: given a set {5\{L)} for some L, the most likely value 
of X is the one that maximizes the likelihood C = Y\iPL{^l{L),x). In practice we work with 
log C We tested this method on an ensemble of 300 spectra, with = 1000 after depletion, 
with values of x = 0.00, 0.01, 0.02, . . . , 0.29, 0.30. The mean, x and standard deviation, ax, 
of the 300 values of x were returned for L. In Fig. H] we see some representative examples 
of log£, and in Fig. [5] we see the results of x for 8 spectra, 4 with x= 0.04, and 4 with 
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x=0.00 




FIG. 3: (color online) The effect of depletion on p{5) is shown here, where we have plotted the 
ensemble average of p{5). The set of 3 lines on the left is for 6 = 53(10), while 5 = (5|(40) is on the 
right. Each set of 3 lines corresponds to a; = 0.00, 0.04, and 0.08. The parameterized p{5,x), Eq. [7] 
is plotted in black , and is barely distinguishable from the ensemble average, in blue. 



0.12. Each set of 4 spectra were randomly chosen, but they are representative of the general 
behavior of x vs L. In Fig. |6]we see the ensemble average for x=0.04, 0.06, . . . 0.014, 0.016. 
We used as errorbars in this plot. The figure suggests that the most reliable range of L 
to use has 20 < L < 40, because in this range, settles down to a smaller value, see Fig. [TJ 
and X is close to the true value of x. 

The MLM doesn't have an error bar for the most likely value of x it returns for a specific 
spectrum. In the analysis of an individual spectrum, one may report a graph of log C vs 
X, and state its maximum. If the peak in log£ is sharp, one would have more confidence 
in the results. Agvaanluvsan et al 12|, used the broadness of the graph of log£ vs x to 
give a range for x. However, if the spectrum being analyzed is from a known ensemble then 

as described above would be a reasonable gauge of how close to the true value of x the 
MLM gets. Looking at ax averaged over either x or L, we are justified in using a value of 
a = 0.006 in our analysis, see Fig. [71 

Spectrum size is an issue of great practical importance when applying these results. In 
the case of neutron resonance data, the spectra we analyzed had typically 80 to 90 levels, 
and we also looked at subsets of them. The acoustic spectra we examined had ~ 250 levels. 
So regardless of the way we choose the error bar, we need to state the dependance of it. 
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FIG. 4: The results of the MLM calculation for 10 randomly chosen depleted members of the GOE 
with X = 0.00, 0.04, and 0.08. The left hand side has L = 20, and the right hand side has L = 40. 

Each spectrum in our test had N = 1000, and this yielded N ~L — 1 values of 6. We suggest, 
based on the behavior of ax that we just calculate x for 10 < L < 40, take ax ~ 0.006, and 
for spectra of size A^, use a = 0.006 y^^- This is obviously a very rough rule of thumb, 
don't forget that when = 80 the lowest values x can have are 1.25%, 2.5%, 3.75%, and 
5%, corresponding to 1, 2, 3 and 4 levels missed. We will see that the drift in the returned 
value of X for a given spectrum is often the biggest consideration for extracting x. 



V. THE BOHIGAS EXPRESSION FOR A3(L,x) 

In this section we use Eq. IHto extract x from the depleted GOE, and compare the results 
with that of our MLM. A comparison of A3(L) for the depleted GOE with Eq. H] gives an 
excellent agreement. In Fig. [8]we see the results for x = 0.00, 0.05, 0.10, 0.15 and 0.20. The 
Bohigas result lies very close to the GOE results in blue. To test Eq. |4]as a tool for gauging 
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L 

FIG. 5: The results of the MLM calculation for randomly chosen spectra with x = 0.04 and 0.12. 

X we get A3(L) for a depleted spectra, and find the x that minimizes Y^ii^^i.^) ~ ^^.{^1 
Repeating this for 1000 spectra for with x from 0.00 to 0.15 in steps of 0.01, we got the 
results are shown in Fig. [9l The mean value of x was very accurate. The error bars, however, 
are much bigger than for the MLM. In Fig. [7] top panel the dashed line is the from this 
analysis. It is much bigger than ax for our MLM. In Table U we include a column of results 
from this method. 

VI. APPLICATION TO NEUTRON RESONANCES AND ULTRASONIC SPEC- 
TRA. 

In Ref.Q,0] the A3(L) statistic and the MLM of Agvaanuvlusaan et al. was used to 
to gauge the completeness of neutron resonance data. The results from both methods were 
consistent with each other. The uncertainties in x for each of these methods was around 
0.03. Here we do a new analysis of some of the same datasets with the new mlm, and report 
the results. A summary of the data sets is in Table [B We analyzed neutron resonance 
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FIG. 6: The results of the MLM calculation for x = 0.02, 0.04, . . . 0.18 and 0.20, for ah L. cj^ was 
used for errorbars. We used this graph as a guide for using ax = 0.006 for 20 < L < 40 for all x. 



data from 7 isotopes in all. The data was taken from the Los Alamos National Laboratory 



website 



24j . Some of the spectra examined did not yield a fiat x vs L graph. In these cases 



the average x is meaningless, and we report that the method is inconclusive. When a; vs L 
is fiat, we report the result of the MLM as x = x ± cr^;, where the average is taken over the 



range 20 < L < 40, and = O.OOGylOOO/D — 30, where D is the number of levels. 

The cumulative level number gives the first indication of the purity of the data. Kinks in 
N{E) leading to smaller slopes would suggest a section of data where levels were missing. 
Sometimes data sets were compiled from different laboratories. The nuclear level density is 
essentially constant in the range of energies of neutron resonance data, so abrupt drops in 
the level density suggest experimental issues. Using this as a guideline, some data sets were 
split into subsets. 
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FIG. 7: The mean value of a averaged over L, top panel, and x, bottom panel. The dashed line in 
th eiupper panel is the result for the analysis based on Eq. H] 

A. 158Gd 

In panel a) of Fig. [10] we see the results of the MLM. There were 93 levels in all, and 
a drop in p{E) at the 60^^ level indicated that the lowest 60 levels were a more pure set 
than the higher 32 levels. The top dashed line, with x = 0.30 is for the higher 32 levels. 
In our MLM, the maximum value we went to was x = 0.30. A spectrum with x > 0.30 
would ideally return a flat line for x vs L at x = 0.30, suggesting in this case, that there 
were originally at least 47 levels in this range. The results for the full set (solid line) are 
X = 0.121 ± 0.024 which translates into there being 93/(1 — x) = 105 ± 3 levels initially, 
13 ± 3 being missed. This is consistent with the lower 60 levels being pure, and the upper 
32 levels being 0.78 of the full spectrum in that range. 
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FIG. 8: A comparison between A3(L) vs L for the Bohigas expression and the ensemble averages. 

B. ^^Ni 

Guided by M{E) we took the the fuU set of 63 levels for the ^'^Ni data. A plot of x vs L 
shown in Fig. [TO] panel b), raises serious questions about the reliability of the MLM in this 
case. The results are inconclusive. 

C. 152gjn 

An analysis of the full set of 91 levels gives an x = 0.081 ± 0.024, as seen in shown in 
Fig. [To] panel c), solid line. This corresponds to there being 99 ±3 levels in the full spectrum, 
with 8 ± 3 missed. The first 70 levels look pure, so it seems that the missed levels were in 
the upper range. 
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FIG. 9: The test results for using the Bohigas expression to extract x. 



D. 234u 

The 118 levels in the full set of resonances yielded higher values of x as L increased, 
solid line in panel d), Fig. [101 As in the case of ^^Fe and ^^Ni, any conclusions are therefore 
suspect. There was a kink in ^{E) after the 78^^ level. These first 78 levels had a monotonic 
increasing x vs L curve, while the top 20 levels had a decreasing curve. Little can be 
concluded from this. 

E. 236u 

There were 81 levels in the ^36^ gg^_ ^jj^g f^jj gg^- j^g^^ _ q_]^24 ± 0.027 corresponding 
to 10 ± 2 levels missed. Guided by a kink in N'{E) at level 70, we analyzed the lowest 69 
levels and got x = 0.038 ± 0.030 corresponding to 3 ± 2 levels missed in that range. 
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F. 235-y 



This is a spectacular data set, with over 3100 levels. Guided by a level density plot, we 
analyzed the lowest 950 levels. The target is odd, with j = |, so the neutron resonances 
were compound states of the ^ss-q nucleus with j = 3 and 4. These resonances were labeled 
with angular momentum, and we separated the 2 sequences of levels and analyzed them 
separately. The result is in panel f). Fig. [TOl where the j = 3 subset is the solid line, with 
X = 0.053 ± 0.010, corresponding to 20 ± 4 levels missed. The dashed line is the j = 4: 
subset. The mean value of x is 0.029 for the j = 3 set, and 0.031 for the j = 4: set, with 
X = 0.045 ± 0.008, corresponding to 26 ± 5 levels missed. Note the ran ge o f L, and how flat 
the lines are. This result is consistent with the other estimates of x in 15|. 



G. Acoustic data 



Here the spectra were resonant frequencies of aluminum cavities. The full experiment is 
described in [18!]. The cavities used in the experiment are made out of aluminum cubes with 
a cube side size of d=20mm. The symmetry of the cube is broken by additional features such 
as an asymmetrically placed cylindrical well and a removed side corner. The radius of the 
well is 5mm and its depth is 18mm. Different configurations of the transducers on the blocks 
gave different data sets. The experimentalists used a comparison of A3(L) with depleted 
GOE results, and conclude that there were 25% of resonances missed. Only data sets 1,2 
and 4 of the 6 data sets analyzed gave x vs L curves less than 0.30. These are plotted in 
Fig. [m The MLM result for all the data sets from 1 to 6 respectively are x = 0.204 ±0.012, 
X = 0.221 ± 0.014, X > 0.30, x = 0.234 ± 0.014, x > 0.30, and x > 0.30. The corresponding 
results for the Bohigas method are 0.28, 0.20, 0.4, 0.23, 0.29 and 0.4 



VII. INTRUDER LEVELS 



In spectroscopy it is quite possible to falsely label a background noise peak as a level, 
or to have a level with a different angular momentum to appear. It is important in RMT 
to know how many sequences of levels are present, (a sequence is a set of energies with the 
same quantum numbers). If an intruder from a different sequence is present, it will not be 
repelled by the other levels, and cause differences in the fluctuation properties of the spectra. 
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FIG. 10: The results of the MLM calculation for neutron resonance data. In panel a) we have 
^^^Gd. The results for the full set (solid line) lies between those for the subset containing just the 
lowest 60 levels, (lower dashed line) and the subset containing the top 32 levels (upper line). In 
panel b) we have ^^Ni, nothing can be concluded from this about x. In c) we see the full spectrum 
of ^^^Sm (solid line) above the lowest 70 levels (lower dashed line. In d) wc have the full ^^^U data 
set (solid line) looking more incomplete that the first 75 levels alone, (middle line dashed), while 
the levels 76 to 118 (lower dashed line) look like a complete subset. In e) we have the full ^^^U set 
of 81 levels looking more incomplete than the lowest 69 levels. Finally in panel f) we see the ^^^U 
data. The j = 3 subset is the solid line, and the dashed line is the j = 4 subset. Note the range 
on the vertical axis. 

In neutron resonance data, s-wave neutrons on an even- A target give a set of spin-^ levels, 
one could have a spin-| level in their midst from a p-wave neutron. An examination of 
jV(5) was performed for spectra with a fraction y of intruder levels. An ensemble of 1200 
spectra of size (1 — y)D were prepared, with D=1000, and were then polluted by adding yD 
"levels" , which were just random numbers with a uniform distribution over the range of the 
spectrum. The calculational details were much the same as for depletion. The results were 
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FIG. 11: The results of the MLM calculation for acoustic data for sets 1,2, and 4. The other sets 
had X > 0.3. The mean values for the interval 20 < L < 40 are 0.20, 0.22 and 0.23 respectively. 

surprisingly similar to those for depletion. In Fig. [121 lower panel, we see p{5) for the cases 
of X = 0.08 depletion, and y = 0.08 intruders. In both cases L = 40. The upper panel shows 
the parameters in the fit to Eq. [5] for both cases. 



VIII. SHELL MODEL SPECTRA 



Recent developments have called into question t 
the fluctuation properties of complex spectra. In 



le vahdity of applying RMT to describe 



17( 1 deviations from the Porter-Thomas 



distribution for reduced neutron widths of s-wave resonances was revealed. This issue was 



addressed in 19[ where the energy dependance of the widths near a maximum of the neutron 
strength function was found to differ from a/E. In 20[ deviations from the FT distribution 
were seen to lead naturally from the a careful description of unstable quantum states with 
open decay channels. The microscopic physics of reactions not captured in RMT was shown 
directly to lead to deviations from the FT distributions [21I where the continuum shell 
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TABLE I: The results for x, the percent of missing levels in the data. 



Isotope 


NND 


A3(L)(Bohigas) 


Pi5) 


N (# levels) 


subset 


58^ 


0% 


18% 


Inconclusive 


63 


All 




3% 


0% 


0% 


91 


1 ^ 70 




3% 


10% 


8% ± 2% 


91 


All 




11% 


13% 


12% ± 2% 


93 


All 




0% 


0% 


0% 


93 


1 ^ 60 


158Gd 


12% 


42% 


>30% 


93 


61 ^ 93 


234U 


9% 


40% 


Inconclusive 


118 


All 


234U 


6% 


13% 


Inconclusive 


118 


1 75 


234 u 


7% 


4% 


Inconclusive 


118 


76 118 


236^ 


5% 


20% 


12% ± 3% 


81 


All 


236 u 


0% 


5% 


4% ± 3% 


81 


1 ^ 69 


235U j = 3 


3% 


9% 


5% ± 1% 


1436 


1 ^ 381 


235u j = 4 


2% 


4% 


5% ± 1% 


1732 


1 569 



model 
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23| was employed. It is reasonable to see if the A^lL) statistic can discriminate 



between the GOE and a model which includes more physics. The shell model with only 
2-body interactions fits the bill. It allows us to get large pure spectra, and include physical 
restrictions. The following calculations were carried out with 12 valence in the sd model 
space with the "USD" interaction of B. H. Wildenthal using the Oxbash code. There are 
5768 levels with J+T = 2+1 and 3276 levels with J+T = 2+0 (see [8:] for more details). 
In FigJTS] we see the well established ^] result for A3(L) for the shell model, it is well 
described by RMT. The N'iS) is well within the bounds set by the variance from one spectra 
to another, and agrees well with the RMT. 



IX. CONCLUSION 



A maximum likelihood method was devised to gauge the incompleteness of experimental 
spectra when a RMT analysis is appropriate. The method is based on the definition of the 
A3(L) statistic. The distribution of random numbers {6}, the mean of which is A3(L), was 
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FIG. 12: Lower panel p{5) for spectra with x=0.08 depletion (left peak) and y=0.08 intruders. 
They are quite close together and it is hard to distinguish the effect of missed levels from intruders. 
The ensemble average data and the fits are superimposed, and indistinguishable. The upper panel 
shows the parameters a, b, and c for in the fit Eq. [71 The solid lines are for intruders and the 
dashed lines are for missed levels. 

parameterized. The cumulative distribution M{6) was accurately fitted with a simple three 
parameter function of log 5: Af{6) = |(1 — Erf[a + 6 log 5 + c]. These parameters, a, b, and 
c were parameterized as functions of x, for each L, yielding a probability density p{6, x) for 
5 with continuous parameter. Our MLM is based on this p{6,x). 

The method was tested on a depleted GOE and returned accurate values of x. Exper- 
imental data was then analyzed. The results for some neutron resonance data sets was 
consistent with earlier analysis, but occasionally no conclusions could be drawn about the 
completeness of the data. The acoustic spectra of an aluminum block was then analyzed. 
The results in 3 out of 6 samples were inconclusive, and for the remaining three, the results 
were consistent with the conclusions on the experimentalists. 

The expression EqJUwas tested and used to gauge x. It was found to give a good ensemble 
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FIG. 13: The A3(L) and M{6) for shell model spectra with J+T = 2+1 (dashed hne), and J+T = 
2+0 (dotted hne), compared to the GOE resuh (sohd hne) for pure GOE spectra. 

average for x but the spread ax was large. The neutron resonance data and the acoustic 
data were analyzed with this expression, and agreed most of the time with the MLM results. 

The question of intruder levels was addressed and the effects on the A3(L) statistic, as 
well as p{6, x) were seen to be very similar. 

The shell model provides us with a long pure sequence of pure levels from a system 
governed by a hamiltonian distinctly different and more physical than those of RMT. Nev- 
ertheless, shell model spectra are well described by RMT [sl . We see that the A/'(5) was no 
exception, and couldn't distinguish the shell model spectra from the GOE. 
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